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Abstract 

In hot plasmas, a temperature of a few tens of eV is sufficient for 
producing highly stripped ions where multipole transitions become im- 
portant. At low density, the transitions from tightly bound inner shells 
lead to electric-quadrupole lines which are comparable in strength with 
electric-dipole ones. In the present work, we propose analytical formulas 
for the estimation of the number of E2 lines in a transition array. Such ex- 
pressions rely on statistical descriptions of electron states and J-levels. A 
generalized "J-file" sum rule for E2 lines and the strength-weighted shift 
and variance of the line energies of a transition array n£ N+1 — > n£ N n'£' 
of inter-configuration E2 lines are also presented. 

I Introduction 

In atomic spectroscopy, the transition probabilities of electric-quadrupole (E2) 
lines are low, as compared to electric-dipole (El) ones. Indeed, the ratio be- 
tween E2 and El transition probabilities is proportional to (ao/XZ*) 2 , where ao 
is the Bohr radius, Z* the effective charge of the plasma and A the wavelength 
PQ . This quantity is very small for most atomic transitions located in the range 
of UV to IR. However, they are important in astrophysics; for instance, Charro 
et al. performed a careful study of intensities of E2 spectral lines in the ion Mg 

II along the spectral series [2 13] . In astrophysical plasmas of Nebulae and of 
the interstellar medium, the density n e of the free electrons, which are respon- 
sible for collisional excitations and de-excitations, is generally very low, and the 
mean time between collisions is longer than radiative E2 lifetimes, even though 
these are much longer than for El transitions. Consequently, whenever an atom 
gets excited by a collision onto a level, metastable or not, its only decay chan- 
nel is by radiative transitions (forbidden or allowed), and forbidden transitions 
appear with an intensity comparable to the allowed ones [1] . Fournier et al. [5 
identified by ab initio calculations a bright 3c? 8 — > 3d 7 4s E2 line for Mo XVII. 
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They also studied the ratio of 3d — 4s E2 to 3d — 4p El emission in Mo XVI as a 
function of density, and noticed that the sharp falloff in the values of this ratio 
for Mo XVI above certain densities provides a plasma diagnostic for planned 
fusion reactors (high-Z impurities breed the dilution of the fusion "fuel" in reac- 
tors and can have a strong impact on spatial current distributions and radiative 
patterns in the plasma) [SJ [7] . Atomic data of tungsten are strongly needed for 
identification of emission lines in future fusion reactors (ITER) where it is used 
as a plasma- facing material [5] . Then, forbidden lines are of great interest for 
plasma diagnostics because the associated radiation intensity is very sensitive 
to density and temperature. Quinet gave a theoretical survey of forbidden tran- 
sitions (E2 and Ml) in the Ap k and id k ground configurations of ions Rh-like 
W 29+ to Ga-like VF 43+ [9] . Recently, Neu et al. pointed out a very strong E2 
spectral line in the soft X-ray region at a wavelength of 0.793 nm originating 
from Ni-like tungsten VF 46+ [TO]. Clementson et al. measured this E2 line to- 
gether with M3 (magnetic octupole) ground state transitions in the same ion 
using high-resolution crystal spectroscopy at the electron-beam ion trap (EBIT) 
facility in Livermore [TT|. EBIT light source was also used by Ralchenko et al. 
[H HI to record spectra of Br-like W 39+ to Co-like W A1+ in the 12-20 nm 
region. On the other hand, in most laboratory plasmas, n e is noticeably higher. 
For medium-range n e , the probability for collisional depopulation may be sur- 
passed only by the probability of radiative decay through allowed transitions. In 
this case, quadrupole transitions are not observed, justifying their name "forbid- 
den". However, they can be important in some circumstances [U] [T5] [TB] . For 
instance, electric-quadrupole decays were observed in the spectra of neon-like 
[TT] and nickel-like [T5] ions produced by laser irradiation. 

As mentioned before, forbidden lines can arise from metastable levels of ex- 
cited configurations 19 . In highly ionized atoms having ground configurations 
3p m the lowest excited configuration is 3p m_1 3d. For m > 1, p m ~ 1 d has levels 
with J' greater by two than the largest J of p m , and downward electric-dipole 
transitions from those levels are therefore forbidden by the dipole selection rule 
|AJ| < 1. Magnetic-dipole and electric-quadrupole transitions within the con- 
figuration Fe IX 3p 5 3d and arising from such metastable levels are responsible 
for several lines of the solar corona. It is interesting to mention that in Pb I, 
parity-forbidden transitions for electric-dipole radiation 6p 2 — > 6p7p have been 
observed [T9J. In that case, the upper levels are not metastable, because El 
transitions to 6p7s are possible. In absorption, electric-quadrupole transitions 
of the type s — > d have been extensively studied in the alkalis (20] [21] . 

The effective core charge (or screened nuclear charge) seen by an electron 
is the net charge Z of the nucleus with the N c — 1 electrons of the ion core 
Z c = Z — N c + 1. Along an iso-electronic sequence, the E2 radial matrix 
elements decrease like Z~ 2 , rather than like Z~ x as do the El matrix elements. 
For inter-configuration transitions with An ^ 0, one has A~e 2 /Ae 1 oc Z 2 , where 
and A^ 2 represent the transition probability rates for El and E2 lines 
respectively jTHJ [23 - For inter-configuration transitions with An = 0, the orders 
of magnitudes of the El and E2 transition probabilities vary in the same way 
with respect to the effective core charge Z c [19] . 

The detailed calculation of all the line energies and strengths in complex 
atomic spectra is difficult and, in some circumstances, useless. Indeed, when 
the densiy is sufficiently high so that the physical broadening mechanisms (e.g. 
Stark shifts) are important and/or when the number of lines becomes large, 
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the lines coalesce into broad structures. Statistical methods (23J 1211 
I2T1 1251 [251 1501 I5T1 1521 155] are required because, experimentally, some quantities 
cannot be determined individually, but only as weighted average quantities (for 
instance, in emission/ absorption spectra of highly ionized atoms). Moreover, 
explicit quantum calculations can be unappropriate, e.g. if the hamiltonian 
matrix is huge. In addition, global methods can reveal physical properties hid- 
den in a detailed treatment of levels and lines. For example, they have opened 
the way to the definition of the generalized J-file sum rule [53] and of effective 
temperatures for the configuration populations |35] . 

Even for El lines, there exists no general compact formula for the exact 
number of lines in a transition array. Using group-theoretical methods, Krasnitz 
[361 137j obtained such a compact formula only in the simple case of configura- 
tions built with non-equivalent electrons. The statistics of El lines was studied 
by Bauche and Bauche-Arnoult [251 EFJ) but very few results were obtained 
as concerns the statistics of E2 lines in complex atomic spectra. The purpose 
of the present work is to provide an analytical formula for the estimation of 
the number of E2 lines in a transition array. Such a quantity is important for 
opacity codes, for instance, in order to decide whether a transition array can be 
described statistically [2H [251 123 123 or requires a DLA (Detailed-Line Account- 
ing) calculation, relying on the diagonalization of the Hamiltonian [315]. Like for 
the El statistical description is used, because this problem does not lend 

itself readily to exact calculations. The first step consists in determining the 
statistics of the angular quantum number J. However, due to the fact that the 
quantum number J is the eigenvalue of no simple operator, its mathematical 
study is tedious. Therefore, it is more appropriate to study the distribution of 
M, the eigenvalue (in units of K) of the operator J z (projection along the z— axis 
of operator J). The J values can be obtained from the M values by means of 
the method of Condon and Shortlcy 1401 . Any unknown distribution can be 
characterized by its moments. The n th moment of the M distribution reads: 

M« = - E (lJM\J z hJMY\ (1) 
q L — ' 

where the sum runs over all the g states of the configuration in intermediate 
coupling. The distribution P(M) of M being symmetrical, the odd-order mo- 
ments are zero. In general, it is admitted |41| that the first four moments are 
sufficient to capture the shape of a distribution. 

In section [2J the second and fourth-order moments of P(M) for a non- 
relativistic configuration are recalled. Section [3] contains the calculation of the 
number of E2 lines between two configurations in the case of three different 
modelings of P(M): Gaussian, Gram-Charlier and Generalized Gaussian. Sec- 
tion [4] contains the calculation, using the aforementioned three distributions, of 
the number of E2 lines inside a configuration. In section [5J a generalized J-file 
sum rule for E2 lines is given, and in section [5] the expressions of the strength- 
weighted shift and variance of a transition array of inter-configuration E2 lines 
are provided for n£ N+1 — > n£ N n'£'. Section [7J is the conclusion. 
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2 Variance and kurtosis of the projection of the 
J angular momentum 

The variance v of the states of configuration £ is given by [3U |3H] : 
where 

4l 2 +4l + 3 

U W = j- • (3) 

The fourth-order moment of configuration £ N reads 

M4 (l N ) = N(U + 2- N)[x(£)N(4£ + 2-N)+ y(l)}, (4) 



where 



2£ — 1 

X(i) = 240(161* -1) (40£3 + 84 ^ + 51) (5) 

and 

2/ + 1 

y W = 60(16^ -1) ( ' 16 ^ _ 24£3 ~ U2 + 2U + 9) - (6) 

For a configuration with w open sub-shells £^£^ 2 £^ 3 ■ ■ , one has for the 
total variance: 

w 

and for the total kurtosis: 



i=l 



M4 («^...e) = E^^f) 

if 

+6 E I* (8) 

where the first sum runs over all subshells i of interest, and the second sum 
runs over all pairs of subshells. In the following section, we will use the 
fourth-order reduced moment: 

«4 = -4, (9) 

v 

which can also be put in the form |42) : 

1 w 

"4 = 3 (a 4 + qih) (10) 

w z * — ' 

t=i 

where 
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_ (24 + l)(16lf + 241? + 8lf - 244 - 9) 

60(4^-1) ' 1 ' 

-16£f ~ 16^ + 88£ 2 + 1364 + 3 

(12) 



120(44 - 1) 



and 



^(44 + 2-AQ 

91 44 + 1 ■ ( j 

3 Number of E2 lines between two different 
configurations 

The initial and final configurations are denoted C and C respectively. The 
selection rules for E2 transitions in intermediate coupling are A J = 0,±1,±2 
with J + J' > 2, and A£ = 0, ±2 with I + 4 > 2 (transitions s ->• s' are strictly 
forbidden for isolated atoms, but have been observed by means of Stark-effect- 
induced mixing of n's' with n"p states [13]). Therefore, the number Le 2 of E2 
lines between two different configurations is given by 



L E2 {C-C) = ]T Qc(J)[Qc>(J) + Qc>(J + l) 

J=Jmin 

+Qc>{J - 1) + + 2) + Q C ,(J - 2)] 

+e'(J rain ), (14) 

where Qc(J) (respectively Qc(J)) is the number of the levels of configuration 
C (respectively C) with total angular momentum J and e(J m i n ) is a small 
correction, J m j n being the smallest value of J. In many cases, J m i n = or 1/2 
according to the parity of the number of electrons. For configuration sd, J m in=l 
and for configurations involving orbitals with a large angular momentum J, the 
value of J m in can be larger. A continuous representation of Eq. (|14l) yields 



L B2 {C-C) w / Qc(J)[<3c(J) + Qc'(J + l) 

+Qc (J-l) + Qc(J + 2) + Qa>(J- 2)] d J 

+e'(J mi „), (15) 

where e'(J m i n ) is a border correction. In fact, e'(0) = e'(|) = are values 
compatible with the accuracy of the replacement of a discrete sum by an integral 
(see Appendix A). Using the second-order Taylor development 



Q{J + k) 



k 2 d z Q 
2 dJ 2 



(16) 



one obtains 
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L E2 {C-C) « 5 / Q C (J) 

J -1/2 

The method of Condon and Shortley [I0J enables one to express Q(J) as 

M=J+1 

Q( J )= E (-l) J - M i , (M) = P(J)-P(J + l), (18) 

where P represents the distribution of the angular-momentum projection M. 
For a configuration i^ 1 ^ 2 ^ 3 ■ ■ ■ i^™ , P(M) is determined through the relation 

P NuN2 ,..{M) = (P Nl <g> Pv 2 <g> Pv 3 <g> • • • ® P W „)(M), (19) 
where the distributions are convolved two at a time, which means that 

(P Nt (g>P Nj )(M)= P N,( M ') x P N .(M - M'). (20) 

M'=— oo 

The total number Le 2 {C — C) of E2 lines is invariant under the following 
transformations: 



Qc>(J) + 



d 2 Q c 



dJ 2 



dJ. 



(17) 



(i) The configurations C and C can be interchanged. 

(ii) In C and C , one can replace simultaneously all the subshells by their 
complementary subshells, e.g.: £ by £ 4 ^+ 2 ~ . Therefore, the C — > C array, 
denoted in general £? 1+1 £% 2 -> ^^ 2+1 has the same number of E2 lines as 
its complementary array: 

^e 1 + l-N 1 ^4£ 2 +2-N 2 _^ ^« 1 +2-Af 1 ^4£ 2 + l-Ar 2 ^1) 



(iii) The C — > C" transition array has also the same number of E2 lines as 
its two following semi- complementary array 

£«i+l-JV 1 ^JV 2 +l _^ |44+2-iVi^iV2 / 2 2) 

and 



The first semi-complementary array is deduced from C — > C" as follows: replace 
the l\ subshells by their complcmcntaries, and exchange the £2 subshells. For 
the second one, replace the £2 subshells by their complementaries, and exchange 
the £\ subshells. Complementarity ensures that the variances are unchanged, 
but it is not the case for semi-complementarity. However, since 
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semi-complementarity can be ensured by replacing vc and Vc> by their half sum 
(vc + Vc)/2 |38| . However, the kurtosis does not follow a relation similar to 
Eq. (j2~4"l) . It can only be ensured that the final formula for Le 2 {C — C) be 
symmetrical in C and C", by setting 04 = («4,c + a 4,C")/2- Following [55] . we 
make the assumption that 



dp 



J+1/2 



(25) 



Transition array 


Exact 


Gaussian 


Gram-Charlier 


GG 


d 4 d 3 s 


887 


953 


885 


839 






(+7.44 %) 


(-0.23 %) 


(-5.41 %) 




1015 


1215 


1092 


1045 






(+19.70 %) 


(+7.56 %) 


(+2.96 %) 




110 


124 


116 


80 






(+12.73 %) 


(+5.45 %) 


(-27.27 %) 


P 4 -> P :i f 


110 


124 


116 


81 






(+12.73 %) 


(+5.45 %) 


(-26.36 %) 


-> d 4 ff 


8299 


9550 


8697 


8415 






(+15.07 %) 


(+4.80 %) 


(+0.14 %) 


dV -> d 4 5 4 


296780266 


299065925 


298981701 


293198984 






(+0.77 %) 


(+0.74 %) 


(-1.21 %) 


dS 1 dV* 1 


762544 


897641 


783942 


774216 






(+17.72 %) 


(+2.81 %) 


(+1.54 %) 



Table 1: Number of E2 lines for different inter-configuration transitions cal- 
culated using Gaussian, fourth-order Gram-Charlier and Generalized-Gaussian 
modelings of P(M) and compared to the exact values. An orbital i corresponds 
to £=6. 



Considering the Gaussian modeling of P(M) 

9 1 

; exp 

/2TTV V 2» 



P(M) = -^= exp 1-^1, (26) 



one obtains 



which yields 



Q(J) = — |=(2J + l)exp (-^^- ) ■ (27) 



where <?c represents the degeneracy (total number of states) of configuration C; 
for instance, the degeneracy of a configuration i^i^if 3 ■ ■ ■ reads 

<») 
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It was also suggested in Ref. 38 to use a fourth-order Gram-Charlier modeling 
of P(M): 



P(M) 



9 



'2irv 

In that case, one finds 



exp 



M 2 
~2T 



a 4 



— (3 - 6— 

24 V v v 



Af 4 

2~ 



(30) 



Q(J) 



which yields 
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(2J+1) 



v\/8ttv 

(2J+1) 3 



-10- 



4 w 



(q 4 - 3) 
24 

(2J+ l) 5 
16 v 2 



15(2 J + 1) 

(2J+1) 2 



exp 



8 u 



(31) 



= 16384^^/2 flO> (-6747 + 2018 t; 
+5 a 4 (938 - 124 v + 21 04 (-11 + 2«)))] 
Finally, considering the Generalized-Gaussian modeling of P(M): 



p(m) = 4= " 



2AT(l + i) 



with A = 



r(f)' 



(32) 



(33) 



where v is a positive real number, and T(x) is the ordinary gamma function, 
one has 



Q(J) 



gcv 



/ 2J+1 V 
2A 2 i; r (1 + i) V 2AV^ J 

and the resulting number of E2 lines is 



exp 



2J+1 



2A^ 



(34) 



L E2 (C-C) 



5 gc gcv 3 2~ 

Q4\3 v 3/2 p (I) 

1 



4"(i/- l)(2i/- 1) 
A 2 ^ 



2-^ 

v , 



AT [2 



(35) 



As can be seen in table [TJ the fourth-order Gram-Charlier expansion series 
provides a better agreement in most of the cases, except for d i 1 — > (Pg 1 ! 1 and 
d 5 —> d A g 1 . This is due to the fact that, when a high-<? electron is involved, the 
shape of the distribution P(M) is not quasi- Gaussian anymore, but exhibits a 
plateau which cannot be modeled by Gram-Charlier expansion. On the 
other hand, the Generalized Gaussian can depict such situations (the "door 
function" corresponds to an exponent equal to v = 9/8). Furthermore, the 
second-order Taylor-series expansion (see Eq. ([16])) is not as precise as in the 
El case, since the selection rules here imply that the maximum value of k is 2 
(against 1 in the El case). 
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4 Number of E2 lines in the intra-configuration 

case (C = C) 



( ,C\T\ T\ 0"\'\VP\ TlOTl 






( . i ' ; 1 1 v i - ( ^firliPT - 


GG 


d b 


358 


380 


352 


335 






(+6.15 7o) 


(-1.68 %) 


(-6.42 7o) 


<rp 


15383 


16499 


15446 


14967 






(+7.25 7o) 


(+0.41 7oj 


(-2.70 To) 


p z f 


353 


402 


358 


322 






(+13.88 7oj 


(+1.42 To) 


(-8.78 70 J 


p 6 f 


611 


696 


622 


567 






(+13 91 %) 


f + 1 80 %) 

( [ X . O W /(J 1 


(-7.20 %) 


9* 


903622 


968412 


909392 


892418 






(+7.17 %) 


(+0.64 %) 


(-1.24 %) 


fV 


5437574 


5839649 


5481760 


5373233 






(+7.39 %) 


(+0.81 %) 


(-1.18 %) 




40832855 


43424047 


41104770 


40379223 






(+6.35 %) 


(+0.67 %) 


(-1.11 %) 


/V 


97615 


107207 


98683 


96137 






(+9.83 %) 


(+1.09 %) 


(-1.51 %) 




690870 


746580 


696669 


681553 






(+0.81 %) 


(+0.84 %) 


(-1.35 %) 




2376 


2545 


2064 


2365 






(+7.11 %) 


(-13.13 %) 


(-0.46 %) 




616 


655 


529 


581 






(+6.33 %) 


(-14.12 %) 


(-5.68 %) 



Table 2: Number of E2 lines inside different configurations calculated using 
Gaussian, fourth-order Gram-Charlier and Generalized-Gaussian modelings of 
P(M) and compared to the exact values. An orbital I corresponds to £—8. 



The number of E2 lines inside a configuration can be estimated as 

1 



-1/2 



Q{J) 



(Q(J)-l) + Q(J + l) + Q(J + 2) 



dJ. 



(36) 



Using the second-order Taylor development presented in Eq. (|16[) , one finds the 
approximate expression 



Li 



1/2 



Q(J) 



5 Q(J) 



d 2 Q 
dJ 2 



,dQ 
'dJ 



A Gaussian expression of P(M) (see Eq. (|2"6"]l ) leads to 

3c (-8V2 ^ + 5«?c(2 «c-3)) 



320F 



, 5/2 



dJ. 



(37) 



(38) 



and a fourth-order Gram-Charlier expression of P(M) (see Eq. (|3"0"|) ) to 
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-2048 <?c(a 4 ,c + 5)v 2 



c 



32768^2^ v 5 J 2 



5 ^ 9C 575 {~ 6747 + 2018 v c + 5 <xi, c [938 - 124 v c 



32768V2tt v 5 J 2 



+21a4 )C (2«c-H)]}. (39) 
Using a Generalized-Gaussian expression of P(M) (see Eq. one obtains 



5 g 2 c v^l- 



3/2 

-4 r 1 2 



128A34/ 2 r(I) 

1 



4V>-1)(2^-1) r /_3^ 



X 2 vc V ^/ 



4A^r 



(40) 



Here also the fourth-order Gram-Charlier expansion series provides a better 
agreement in most of the cases, except for configurations with high-i? electrons 
(see table [5]). The results are closer to the exact values in the intra-configuration 
case than in the inter-configuration case. This is due to the fact that the number 
of lines in the intra-configuration case is evaluated with the exact variance and 
the exact kurtosis of the distribution of Al (in the inter-configuration case, they 
were obtained as the arithmetic average over the initial and final configurations) . 

It is worth mentioning that global methods were also used for the study of 
statistical properties of Auger amplitudes and rates [U] and the dispersion of 
gyromagnetic ratios in complex spectra [45 . In appendix C, we give, following 
the work of Kyniene et al. |44) . an approximate expression for the number of 
Auger amplitudes with the three different modelings of P(M). 



5 A generalized J-file sum rule for E2 lines 

Using the second-quantization technique [3S], Bauche et al. [33] established a 
linear relationship between the angular coefficient of the exchange Slater integral 
G 1 (n£,n'£ — 1) and the J-file sums of the El line strengths defined by Condon 
and Shortley [JO]- In the case of E2 lines, a J-file sum rule can be also obtained, 
following the same procedure as for El lines. The sum of the strengths of all 
E2 lines starting from a level 7 J of the upper configuration in the £ N £' N +1 — > 

gN+l£/N arra y jg gi ven by; 



s 



A 2 



(l N t' N ' +1 } 7 J - £ N+1 e' N '] = (2J+1) 



(N' + 1) 



{21' + 1) 



(£\\CW\\£>) 



+C(G 2 ; 7 J) 



[I(n£,n'e')f 



with 



/•OO 

I(n£,n'£')= / R ne (r)r 2 R n , e ,(r)dr 
Jo 



(41) 



(42) 
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and 



|C (2) ||0 = (-1)V(2*+1)(2*' + 1) ( 



£ 2 £' 




(43) 



The quantity C(G 2 ;jJ) represents the coefficient of the G 2 (n£,n'£') Slater in- 
tegral in the electrostatic energy of the level jJ. In the present case, one has, 
if \£' -£\ = 2 and £> = max(f, £'): 

and if £ = £': 

These results remain unchanged if passive subshells are added to both configu- 
rations. It is interesting to mention that the total strength of transition array 
C -> C with C = n£ N+1 n'£' N 'n"£" N " ■■■ and C = n£ N n'£' N ' +1 n"£ N " ■ • • 
reads 



where gc represents the degeneracy of configuration C, i.e. 



(46) 



, 4^ + 2 \ A£' + 2\ [41" + 2 . 

■ 9c= U + i J v iv' J v TV" )•••■ (4,) 

6 Shift and variance of a transition array n£ N+1 —, 
n£ N n'£' of E2 lines 

To our knowledge [37] . Harrison and Johnson [35] are the first who introduced 
the term "transition array" for the entire bunch of lines resulting from tran- 
sitions between two configurations. The first experimental spectrum showing 
transition arrays was published by Edlen in 1947 [49 : it concerned transition 
elements in the XUV range and the spectra were generated by low-inductance 
discharge lamps. Since then, such arrays have been observed in a very large 
variety of spectra (see for instance [501 EB [54]). Although, in tokamak 

applications, E2 lines appear usually as isolated (only a few distinct lines are 
visible), it is also possible to estimate the global properties (strength- weighted 
moments) of a transition array of E2 lines. In the present section, we con- 
sider the case of inter-configuration E2 lines for the specific transition array 
n£ N+1 -> n£ N n'£' with £' - £ = 0, ±2 and £ + £' > 2. The moments of the 
distribution of the line energies can be written as: 

kA Ea b [(b\H\b)-(a\H\a)} n \(a\Q\b)\ 2 

Mn = - — (48) 
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where a and b run over all the exact eigenstates of the Hamiltonian in config- 
urations C and C , respectively, and Q is the z component of the quadrupole 
transition operator. The first two moments were calculated by Bauche-Arnoult 
et al. [23] using the second-quantization techniques of Judd [35]. The most 
useful quantities are the mean energy Mi and the variance a 2 = M.2 — M.\ 
related to the spectral width of the UTA. For the mean energy of the array, it 
was found in Ref. 124 that 



Mi = E c > - E C + 6E, 



(49) 



where Ec and Eq> are the average energies of the initial and final configurations 
respectively. There exists a shift between the weighted average energy of a 
transition array and the difference of average energies of the initial and final 
configurations. The shift SE stems from the weighting factor (a|Q|6)| 2 and is 
non zero for transition arrays n£ N+1 —> n£ N n'£'. It was shown in Ref. [23] that 
a 2 can be written as 



^ ] di(k,k ,£,£ , 
k,k' 



•) x ei(n£,n'£',---) 



(50) 



where arc numerical coefficients depending on the number of equivalent elec- 
trons N, di(k, k',£, £',■■■) are combinations of 3nj (n—1,2 and 3) symbols in- 
dependent of N and ei(n£,n'£', ■ ■ ■) are products of Slater integrals of ranks k 

represent the direct Slater integrals in 
the exchange Slater integral of 

C . As in Ref. [23] , we define the quantities: 



and k' . In the following, Fq 



and Fq, 



configurations C and C, respectively, and G c , 



X 


= N(N + 


l)(i£- 


N)(4£-N + l) 


= ( N 


+ l)(4i-N)w, 


V 


= N(N~ 


l)(i£- 


N + 1)(A£ — N 


+ 2) = 


(N- 1)(4£- N + 2)w 


z 


= N(N — 


1)(4£- 


N)(U-N + 1) 


= ( N 


- l){U-N)w, 


u 


= N(4£ - 


N)(M- 


-N + l) = {A£- 


- N)w, 




V 


= N{N — 


1){U- 


N + 1) = (N- 


l)w, 




w 


= N(U - 


N+ 1) 









(51) 

In the following, in agreement with the convention of Ref. [24] . indices n and 
n! are omitted in the Slater integrals. The shift is given by: 

^ = ^ ( %i+i) +1) [E/^(«o+E^w ] (52) 



. fc/0 



with 



and 



t £ k £ \ I £ k £ \ \ £ k £ i — > ■ 

h= l H ) i £' 2 £' f ' ' i > ! 



5 6h ' 2 2(2£+\){U> + 1)) ' (M) 







k 




9k = 


it 
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where Sij represents a Kronecker's symbol. The variance can be written as 



<r 2 = 5>. 

As compared to Ref. [21], the term Hi remains unchanged: 



Hi = 



and one has 



EE 

k=£0 fe'#0 
— ( — l) fc+ ' £ 



V(2fc+1) (2£+l)(4£+l) 

(2* + l) 3 







I 


k 




{ 


£ 


I 


fc' 


}) 


I 




\ 


I 


fc' 





) 











(4^-l)8£(4£+l) 

2 



(55) 



xF*(U)F${U) + yF&(e£)Fg,(U) - 2zF*{U)F$,(U) , (56) 



H 9 



ee(«« 



V(2fc + 1) 



. fc i' £ I k' 



e i 1 

i i 



-i-iy 



iik 



i 



(2* + l) (4€+l) 



ilk' 

i 1 e 2 

iik 1 ' 
e e 2 



(2^+l) 3 (2f + 1) 
(4^-1)2^+1) 





fc 






fc' 






fc' 


















5) 


(J 









x uF£(ll)F£,(ll') + vF£,(ll)F£,{ll') 



(57) 



EE 

fe^O fc' 



1 



fc fc' 2 
I' I I 



6 k 



(2l + l)(M + l) V5 

2 



fc' 2 



£ fc' £' 



I k I 




i i 

(2^ + 1) 

2 





x uF^{ll)G k c, {U') + vF$, (11)0% (11') 



I 1 I 1 k 
Ilk 1 



(2l+l) 3 (2l' + l) 
(4^-1)2^+1) 



(58) 
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Hi 






fc 




o 





i ( 



(2l+l) 2 (2l' + 1) 2 
4^(4^+1) 



£' fc' £' 




xwF£,(££')F$, (££'), (59) 



EE 



2<5fc,fc' 



(2i' + l)(2fe + l) 



(4£+ 1) V5 



r<5fc,2 - 



1 



(21' + 1) 



£ £' k 
£' 2 £ 
fc' £ £' 



2 L 

5 fc '' 2 2£' ^ 



(2l+l) 2 (2f + 1) 2 f £ k £' 
X 4^(4^ +1) ^000 

x«;G^(«')G^(«'), 
and the last term (for the electrostatic part): 







(60) 



Ha 



EE 

fc^O fc' 

-2(-iy 



(-l) fe 




(2f + 1) 




r fc fc' 


2 






* 2), 




£ fc / 1 





fc' 2 



fc fc' 2 

1 

' (2f + 1) 



(4* + 1) 
(2£ + l) 2 (2£' + l) 2 



f fc f 




£ fc' £' 




2£(4^ + 1) 

2 





xwF£,(££')G h c, (££'). 

The spin-orbit contribution is the same as for the El case 



(61) 
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H 7 



-N 



+N{U -N + 2) 



(N + l)(4l- N + 1) 



N(U-N + 1) 



(U-N + l) 



£(£ + !)+£'{£' 




(62) 



where £c (n£) and Cc ( n ^) are the spin-orbit integrals in the respective configu- 
rations C and C. Higher moments can also be calculated [251 I2H1 1301 15T1 155] . 
which requires to use another modeling function than the Gaussian. The Gen- 
eralized Gaussian method [SB] can be applied, but the Normal Inverse Gaussian 
was shown to provide a better depiction of the array when the first four moments 
are known [53 [58] . 

The range of energy of the levels of the upper configuration C responsible 
for the preferential emission constitutes an "emissive zone" [33]. In the same 
way as for an El UTA, the shift and width of the emissive zone of an E2 UTA 
are deduced from those of the complete UTA by restricting them to the radial 
parameters related to the upper configuration of the array. 

7 Conclusion 

We propose analytical expressions for the number of electric-quadrupole (E2) 
lines both in the inter- and intra-configuration cases. The resulting formu- 
las are based on three different modelings of the distribution of the angular- 
momentum projection M (Gaussian, fourth-order Gram-Charlier and Gener- 
alized Gaussian). The Gram-Charlier modeling gives satisfactory results but, 
for high-f electrons, the Generalized-Gaussian distribution is more accurate (as 
for El lines, see Ref. [32). The results are better in the intra- than in the 
inter-configuration case, which can be explained by the fact that the number of 
lines between two different configurations is evaluated with averaged parameters 
(variance and kurtosis of the distribution of the angular-momentum projection 
M). We also provide the expression of the generalized J- file sum rule giving the 
total strength of the lines arising from a given level of the upper configuration, 
together with the formula for the average strength- weighted position and vari- 
ance of E2 lines. The next step will consist in investigating the statistics of the 
amplitudes and strengths of E2 lines. 

8 Appendix A: Boundary effect on the number 
of lines of an array 

At very small values of J, we can linearize the exponential in Eq. (|27l) . i.e. 
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Qc{J) = K c x (2J+1) (63) 

and 



Qc>(J)=K C 'X(2J+l), (64) 
where Kc and Kc> do not depend on J. 

(i) Inter-configuration case: 

In the inter-configuration case, the error committed in considering that 
e'(|) = in Eq. JT5]) is equal to 



^e, inter — 5 / Q c {J)Qc>(J)dJ 

J -1/2 



1/2 

Using the linearized expressions (|63[) and (|64p . one obtains 



(65) 



A eMet « jk#C + #C'). (66) 

6 



(ii) Intra-configuration case: 

In the intra-configuration case, the error committed in considering that 
e'(±) = in Eq. JTS) is equal to 

^e, intra — 3 / Q c (J)Qc(J)dJ (67) 

J -1/2 

Using the linearized expressions (|6"3"|) and one obtains 

A ei i ntra « . (68) 

9 Appendix B: Estimation of J m ^ 

Usually, J m i n is equal to (for integer values) or to 1 /2 (for half- integer values) . 
However, in some circumstances [22] (for instance in the case of configurations 
containing an electron in an orbital with a high angular momentum £), the 
distribution P{M) exhibits a plateau, for which we showed that the Generalized 
Gaussian in a good aproximation 42 . In that case, J mm differs from and 1/2. 
An estimation of J m i n can be obtained through the relation: 

Q(J) > ~. (69) 



Using the expression of Q(J) given in Eq. ([34]) . we find 

/i _ \ r -i l/f i 

Jmi„~A^ -)W —— — A 77 ^ 1 " --, 

\ v J L 1 — v 2 

with 



(70) 
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A-^r(i), (7i) 

and 14 W[x] represents Lambert's function, solution of x e x — y. The function 
W can be expanded as [55 | 150 1 IBT ] 153] : 



OO OO 



W{x) = In* - ln(ln,) + £ £ ^T^S ( 72 ) 



where 



c fcTn = ^-^-S[fe + m,fe + l], (73) 
to! 

S 1 ^, <?] being Stirling number of the first kind [63]|64], also denoted Sp q \ s(p,q) 



Stirling numbers can be obtained by recursion relations [60] . and an 



explicit expression was provided by Karanicoloff [65] . However, since we are 
only interested in an approximate formula, the first two terms In a; — ln(lnx) are 
sufficient. 



10 Appendix C: Number of Auger amplitudes 

Auto-ionization from a state that involves a hole in an inner subshell of the core 
is known as the Auger [66] effect, and the ejected electron is called an Auger 
electron. The term "auto-ionization" (applied to levels produced by excitation 
of loosely bound electrons) was coined by Shenstone [57J . In the non-relativistic 
approximation the amplitude of Auger transitions is equal to the reduced matrix 
element of the Coulomb interaction operator H c : 

(C^J\\H c \\C'j'J'e£jJ) = V2JTT (CjJ\H c \C'~f'J'e£jJ), (74) 

where C is the configuration of an atom, J is the quantum number of total 
angular momentum, 7 denotes all the additional quantum numbers and e is the 
energy of the Auger electron. In intermediale coupling, the number of Auger 
amplitudes for a given channel e£ or egj can be obtained as the number of reduced 
matrix elements of the scalar operator acting between two configurations: 

/>oo 

AWrCC - C) w / Qc{J)Qc»{J)dJ, (75) 
Jo 

where C" denotes Gel and C is the final configuration of the ion. In the 
configuration C" , the contribution of the Auger electron to the variance and 
kurtosis of the distribution P(M) is equal to the one of an electron of the 
discrete spectrum v{e£) — v(i), where v(£) is defined in Eq. ([3]) and 

A* 4 (e^) = 7^( 2 ^ + 1)(48^ 4 + 961 3 + 1521 2 + 104£ + 15). (76) 
It was shown by Kyniene et al. |44] that the Auger transition 

e U 1+ l e N H N 3 ^ ^+2 e N 2 -l e N 3 -l ee (77) 
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does not change under the replacements N2 — > 4^2 + 4 — N2 and N3 — > 4£s + 
4 — N3. In the same way, the transition 

gUi+lg, e U 1+ 2 e N 2 -2 el (7g) 

does not change under the replacement N2 —> 4^2+4 — N2. This is a consequence 
of complementarity, and can be explained by the fact that the recoupling does 
not change the number of matrix elements [33]. Here also, since 



= v (e^ +1 ) + v (£f=+ 4 -^) + « (£f 1+2 ) + v (f A 2 l2+2 - N2 ) + v(e£), 

(79) 

we set v — (vc + t>c ) /2 and, in order to obtain an expression symmetrical in 
C and C, a 4 — (0:4,0 + &4,C')/2- Using the Gaussian expression of P{M) (see 
Eq. (1211), Kyniene et al obtained [33]: 

N Augcl .(C-C')^^0j,. (80) 

The Gram-Charlier modeling of P(M) (see Eq. (J30) leads to 

at tn r>\~ 1009 + 5 a 4 (-62 + 21 q 4 ) 

- C ) « 5C 5C 8192 ^, 3/2 . (81) 

and using a Generalized-Gaussian approximation of P(M), we find 

AW(C ~ C") « 3C ^ 22+1/ , Ar [ 1+ l )0J - (82) 
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